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Abstract 

We present a measurement of the hadronic invariant mass squared (M^) spectrum in charmed 
semileptonic B meson decays B — > X c £u based on 140 fb" 1 of Belle data collected near the 
T(4S) resonance. We determine the first, the second central and the second non-central moments 
of this spectrum for lepton energy thresholds ranging between 0.7 and 1.9 GeV. Full correlations 
between these measurements are evaluated. 

PACS numbers: 12.15.Hh,14.40.Nd,13.25.Hw 
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I. INTRODUCTION 



Inclusive semileptonic decays of B mesons to charmed final states provide an avenue for 
measuring the Cabibbo-Kobayashi-Maskawa (CKM) matrix element |V^>| [l[ and for deter- 
mining non-perturbative hadronic properties of the B meson. In particular, the moments of 
the hadronic mass in B — > X c lv decays calculated in the framework of the Operator Product 
Expansion (OPE) and the Heavy Quark Effective Theory (HQET) 0, H, 0, lEf depend on 
the 6-quark mass (nib) and a few non-perturbative matrix elements that also appear in the 
expression of the total semileptonic width. Thus, measurements of the hadronic invariant 
mass moments 0, Q, @, 0] allow the determination of these non-perturbative parameters from 
the data and reduce the theoretical uncertainty in the extraction of \V c b\ from measurements 
of the semileptonic branching fraction. An improved knowledge of m& also results in a more 
precise determination of \V u b\ from inclusive charmless semileptonic B decays. 

This analysis uses T(AS) — > BB events in which the hadronic decay of one B meson is 
fully reconstructed. The semileptonic decay of the other B is inferred from the presence of 
an identified lepton (electron or muon) amongst the remaining particles in the event. We 
calculate the first two moments of the hadronic invariant mass squared (M\) distribution 13 
directly from the measured spectrum after the effects of finite detector resolution have been 
removed using the Singular Value Decomposition algorithm (Tl| . 

The measurement described in thispaper improves the results previously reported by the 
BaBar and CLEO collaborations 0, [7]. The sensitivity to nib and other non-perturbative 
parameters is increased by lowering the minimum lepton energy threshold to 0.7 GeV. Fi- 
nally, this analysis minimizes the dependence on particular B —>■ X c lv model assumptions 
by calculating the moments directly from the unfolded M\ spectrum. 



II. EXPERIMENTAL PROCEDURE 
A. Data Sample and Event Selection 

The data used in this analysis were taken with the Belle detector [l2| at the KEKB 
asymmetric energy e + e~ collider [l3j |. Belle is a large-solid-angle magnetic spectrometer 
that consists of a three- layer silicon vertex detector, a 50-layer central drift chamber (CDC), 
an array of aerogel threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of- 
flight scintillation counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals (ECL) located inside a super-conducting solenoid coil that provides a 1.5 T magnetic 
field. An iron flux-return located outside of the coil is instrumented to detect K\ mesons 
and to identify muons (KLM). 

The data sample consists of 140 fb _1 taken near the T(AS) resonance, or 152 x 10 6 
BB events. Another 15 fb _1 taken at 60 MeV below the resonance are used to estimate the 
non-BB (continuum) background. The off-resonance data is scaled by the integrated on- to 
off-resonance luminosity ratio corrected for the 1/s dependence of the qq cross-section. 

A generic BB Monte Carlo (MC) sample equivalent to about three times the integrated 



luminosity is used in this analysis. MC-simulated events are generated with EvtGen 14 



and full detector simulation based on GEANT3 [15j is applied. The decays B — > D*£u 



and B — > DCv are generated using an HQET inspired form factor parameterization [16 
The decays B — > D**iv [l7| are simulated according to the Leibovich-Ligeti-Stewart-Wise 
(LLSW) model [l8[ (both relative abundance and form factor shape). The B — > X c iv model 
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D^ixtv decays which are generated using the Goity- Roberts 
The model for the B — > X,Jtv background is a hybrid mixture of exclusive modes 

Light-cone sum 



also includes non-resonant B 
model 



and an inclusive component described by the De Fazio-Neubert model [20 



rule form factors 21, 22| are used for B — > iriu, piv and uoiv. Other exclusive modes are 
simulated according to the ISGW2 model j23|. QED bremsstrahlung in B — > X£u decays is 
included using the PHOTOS package 0. 

Hadronic events are selected based on the charged track multipl icity and the visible energy 
in the calorimeter. The selection is described in detail elsewhere 25 . 



B. Full-reconstruction Tag 

We fully reconstruct the hadronic decay of one B meson (S tag ) using the decay modes 
B+ _> J D(*)° 7 r+, D^°p+, D^°af and B° -> D^~^ + , D^'af Pairs of photons 

satisfying E y > 50 MeV in the laboratory-frame and 118 MeV/c 2 < M(77) < 150 MeV/c 2 
(±3.3<r around the 7r° mass) are combined to form 7r° candidates. K s mesons are recon- 
structed from pairs of oppositely charged tracks with invariant mass within ±30 MeV/c 2 
(±5.1cr) of the nominal K s mass and a decay vertex displaced from the interaction point. 
Candidate p + and p° mesons are reconstructed in the 7r + 7r° and 7r + 7r~ decay modes, requir- 
ing their invariant masses to be within ±150 MeV/c 2 of the nominal p mass. Candidate 
af mesons are obtained by combining a p° candidate with a charged pion and requiring an 
invariant mass between 1.0 and 1.6 GeV/c 2 . D° candidates are searched for in the K~n + , 
K~7T + 7r°, fT~7r + 7r + 7r~ , K s 7i + tt^ and K s ii decay modes. The K"-k + -k + and K^7t + modes are 
used to reconstruct D + mesons. Charmed mesons are selected in a window corresponding to 
±3 times the mass resolution in the respective decay mode. D* + mesons are reconstructed 
by pairing a charmed meson with a low momentum pion, D* + — ► D tt + , D + tt° . The decay 
modes D*° — > D°ir and D*° — > D Qr y are used to search for neutral charmed vector mesons. 

For each B tag candidate, the beam-energy constrained mass M^ c and the energy difference 
AE are calculated, 

M bc = ^J{E heam ) 2 - {p B ) 2 , AE = Eb — £ beam , (1) 

where -Ebeam, Pb and E B are the beam energy, the 3-momentum and the energy of the 
B candidate in the T(4S) frame. In M^ c and AE, the signal peaks at the nominal B mass 
and zero, respectively. We define the signal region by the selections Mb c > 5.27 GeV/c 2 and 
\AE\ < 0.05 GeV. If multiple candidates are found in a single event, the best candidate is 
chosen based on the proximity of AE, Md and AM to their nominal values, where Md is 
the reconstructed D meson mass and AM is the difference between the reconstructed D* 
and D meson masses. Without making any requirement on the decay of the other B meson, 
the number of B + (B°) tags in this region, after subtraction of continuum and combinatorial 
backgrounds, is 61, 365 ± 531 (41, 027 ± 368), Fig. d 



C. Lepton Reconstruction 

Semileptonic decays of the other B meson (-B S i gn ai) are selected by searching for an identi- 
fied charged lepton (electron or muon) within the remaining particles in the event. Electron 
candidates are identified using the ratio of the energy detected in the ECL to the track 
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FIG. 1: Afb c distributions for charged and neutral Stag candidates after requiring \AE\ < 0.05 GeV. 
No constraints are made on the signal side. The points with error bars are on-resonance data 
after subtraction of the scaled off-resonance data. The combinatorial background (cross-hatched 
histogram) is estimated using MC simulation. 



momentum, the ECL shower shape, position matching between track and ECL cluster, the 
energy loss in the CDC and the response of the ACC counters. Muons are identified based on 
their penetration range and transverse scattering in the KLM detector. In the momentum 
region relevant to this analysis, charged leptons are identified with an efficiency of about 



90% and the probability to misidentify a pion as an electron (muon) is 0.25% (1.4%) |27|, [28 

We further require electron (muon) candidates to originate from near the interaction 
vertex, have a laboratory-frame momentum greater than 0.3 GeV/c (0.6 GeV/c) and satisfy 
17° < 9 < 150° (25° < < 145°), where 9 is the polar angle in the laboratory-frame relative 
to the beam direction. If more than one charged lepton candidate is found in the event, 
we only keep the one with the highest momentum in the B rest frame. Electrons from 
photon conversion are vetoed by rejecting the event if the invariant mass of the electron 
candidate and another oppositely charged particle in the event is below 0.04 GeV/c 2 and 
secondary vertex criteria are satisfied. If the charged lepton candidate is consistent with the 
decay J/ip (i.e., the invariant mass of the lepton candidate and another oppositely 

charged lepton in the event is between 3 GeV/c 2 and 3.15 GeV/c 2 ), the event is also rejected. 

In B + tagged events, we require the lepton charge to be consistent with a prompt semilep- 
tonic decay of 5 S ignai- In B° events, we make no requirement on the lepton charge. In 
electron events, we partially recover the effect of bremsstrahlung by searching for a photon 
with laboratory-frame energy E 1 < 1 GeV within a 5° cone around the electron direction at 
the interaction point. If such a photon is found, it is merged with the electron and removed 
from the event. 
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D. Hadronic Mass Reconstruction 

The 4-momentum p x of the hadronic system X recoiling against tv is determined by 
summing the 4-momenta of the remaining charged tracks and unmatched clusters in the 
event. We exclude tracks passing very far away from the interaction point or compatible 
with a multiply reconstructed track generated by a low-momentum particle spiraling in the 
central drift chamber. Unmatched clusters in the barrel region must have an energy greater 
than 50 MeV. Higher thresholds are applied in the endcap regions. 

To improve the resolution in M x , we reject events with a missing mass larger than 
3 GeV 2 / c 4 . Further improvement is obtained by recalculating the 4-momentum of the X sys- 
tem, 

Px = (Pe+-beam + Pe--bcam) ~ PB tag ~ Pi ~ Pv , (2) 

taking the neutrino 4-momentum (E u ,p u ) to be (|p m iss|,Pmiss), where p m i SS is the missing 
3-momentum. Defined as the half width at half maximum, the resolution in M x obtained 
from p' x is about 0.8 GeV 2 /c 4 , compared to 1.4 GeV 2 /c 4 in M\ from p x . 

E. Backgrounds in the Hadronic Mass Spectrum 

We consider the following contributions to the background in the M\ spectrum: non-BB 
(continuum) background, combinatorial background, background from secondary or fake 
leptons and B — > X u iv background. Combinatorial background are true BB events for 
which reconstruction or flavor assignment of the tagged B meson is not correct. 

The shapes of these background components in M x are determined from the MC sim- 
ulation, except for the continuum where off-resonance data is used. The shape of the fake 
muon background is corrected by the ratio of the pion fake rate in the experimental data 
over the same quantity in the MC simulation, as measured using kinematically identified 
pions in Kg — > tt + tt~ decays. We derive the shape of the combinatorial background from 
the generic BB simulation by selecting events in which the reconstruction of B tag does not 
correspond precisely to what was generated in the simulation. 

The continuum background is scaled by the integrated on- to off-resonance luminosity 
ratio, taking into account the cross-section difference. The MC-prediction of the combina- 
torial background is normalized to the data using the side-band region (Mt, c > 5.27 GeV/c 2 
and 0.15 < \AE\ < 0.3 GeV). The normalization of the secondary or fake lepton background 
is found from the data by fitting the electron (muon) momentum distribution p\ |29j in the 
B meson rest frame in the range from 0.3 to 2.4 GeV/c (0.6 to 2.4 GeV/c). The X u lv com- 
ponent is normalized to the number of B + (B°) tags, assuming a branching fraction of 
2.08 x 10- 3 (1.92 x lO- 3 ) for B + X Q J + u (B° -> X~l + v) 

The background in the M\ spectrum is estimated separately in the four sub-samples, 
defined by the charge of _B tag (B + , B°) and the lepton type (electron, muon). 

The purity of the B — > X c tv signal depends on the sub-sample and the lepton energy 
threshold, typical values being around 75%. Table [I] shows the numbers of signal events and 
purities for each combination of Stag charge, lepton type and lepton energy threshold. 
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TABLE I: Number of B — > X c tv signal candidates and signal purity in the four sub-samples, as a 
function of the lepton energy threshold. The yields are quoted with their statistical uncertainty; 
the corresponding signal purity is given in parentheses. 



min 


B+ 


electron 




B + muon 


B° electron 


B° muon 


0.7 


4105 ± 100 (70.5 


%) 


3739 ± 108 (61.5%) 


2491 ± 80 (65.9%) 


2400 ± 86 (60.3%) 


0.9 


3855 ± 


95 (73.2 


%) 


3591 ± 104 (64.8%) 


2353 ± 76 (73.4%) 


2307 ± 83 (67.3%) 


1.1 


3466 ± 


86 (74.9 


%) 


3305 ± 96 (68.3%) 


2098 ±68 (77.1%) 


2120 ± 76 (74.2%) 


1.3 


2894 ± 


72 (75.8 


%) 


2857 ± 84 (70.6%) 


1749 ± 58 (80.4%) 


1800 ± 66 (78.0%) 


1.5 


2195 ± 


56 (74.6 


%) 


2225 ± 66 (72.3%) 


1322 ± 45 (84.2%) 


1388 ± 52 (79.7%) 


1.7 


1384 ± 


38 (77.2%) 


1415 ± 44 (72.4%) 


824 ± 30 (83.7%) 


878 ± 34 (80.7%) 


1.9 


571 ± 


19 (73.8 


%) 


627 ± 22 (74.0%) 


353 ± 15 (84.3%) 


376 ± 17 (76.7%) 



F. Unfolding and Moment Calculation 

We measure the M\ spectrum in 45 bins in the range from to 15 GeV 2 /c 4 (bin width 
0.333 GeV 2 /c 4 ), which is shown in Fig. [21 and unfold the finite detector resolution in this 
distribution using the Singular Value Decomposition (SVD) algorithm [ll[. The unfolded 
M\ spectrum has 15 bins in the range from M|> to about 15 GeV 2 /c 4 . The bin width is 
1 GeV 2 /c 4 , except around the narrow states - D, D*, D\ and D\ - where smaller bin sizes 
are chosen. 

The unfolding is done separately in each sub-sample (B + electron, B + muon, B° electron 
and B° muon). From the unfolded spectrum, we calculate the first moment and its statistical 
uncertainty squared, 

«)^, Am)) .^m^k. (3 ) 

Here, x' is the unfolded spectrum corrected for slightly different bin-to-bin efficiencies and X 
is its covariance matrix, also determined by the SVD algorithm. (Mj^j is the central value 
of the z-th bin of the unfolded spectrum. The second central and non-central moments, 
((Mjr — (M|)) 2 ) and (Mj^) are calculated from the same spectrum, substituting M\ by 
{M\ - {M\)f and M| in Eq. El respectively. 

As the hadron mass moments are not expected to depend on the B meson charge or the 
lepton type [3, 0], we take the average over the four sub-sample results. 

We have tested the entire measurement procedure including event reconstruction, un- 
folding and moment calculation on MC simulated events and no significant bias has been 
observed over the full range of lepton energy thresholds. 



8 





FIG. 2: Measured hadronic mass spectrum for different lepton energy thresholds. The points 
with error bars are the experimental data after subtraction of the continuum background. The 
histograms show the B — > X c lv signal and the different background components, explained in more 
detail in the text. 

III. RESULTS AND SYSTEMATIC UNCERTAINTIES 
A. Results 

Our measurements of (M^), ((M^ — (Mj|)) 2 ) and (M|) for different lepton energy thresh- 
olds are shown in Table [TT] and Fig. [31 The sub-sample results for a given charge of B tSug 
(B + , B°) or lepton type (electron, muon) are compatible within their statistical uncertainty 
only. 
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TABLE II: Measurements of (M|), ((M£ - (M^)) 2 } and (My) for different lepton energy thresh- 
olds. The results in this table are the averages of the four sub-samples, defined by the charge of 
-Btag B°) and the lepton type (electron, muon). The first error is statistical, the second is the 
estimated systematic uncertainty. The different measurements are highly correlated (Tables |VTHX|) . 



EU (GeV) {Ml) (GeV 2 /c 4 ) ((M 2 - (M 2 » 2 ) (GeV 4 /c 8 ) (M|) (GeV 4 /c 8 ) 



0.7 


4.403 ± 0.036 ± 0.052 


1.494 ±0.173 ±0.327 


20.88 ±0.48 ±0.77 


0.9 


4.353 ± 0.032 ± 0.041 


1.229 ±0.138 ±0.244 


20.18 ±0.40 ±0.58 


1.1 


4.293 ± 0.028 ± 0.029 


0.940 ±0.098 ±0.137 


19.37 ±0.33 ±0.36 


1.3 


4.213 ± 0.027 ± 0.024 


0.641 ± 0.071 ± 0.080 


18.40 ±0.29 ±0.26 


1.5 


4.144 ±0.028 ±0.022 


0.515 ±0.061 ±0.064 


17.69 ±0.28 ±0.23 


1.7 


4.056 ± 0.033 ± 0.022 


0.322 ± 0.058 ± 0.040 


16.77 ±0.32 ±0.21 


1.9 


3.996 ± 0.041 ± 0.021 


0.143 ± 0.056 ± 0.038 


16.11 ±0.38 ±0.20 



TABLE III: Breakup of the systematic error on (M|). Refer to the text for details. 



A(M X ) (GeV 2 /c 4 ) 



E* min (GeV) 


0.7 


0.9 


1.1 


1.3 


1.5 


1.7 


1.9 


secondary /fake leptons 


0.033 


0.023 


0.013 


0.007 


0.004 


0.002 


0.000 


combinatorial background 


0.006 


0.004 


0.003 


0.002 


0.002 


0.002 


0.000 


continuum 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


B — » X u lv background 


0.004 


0.004 


0.004 


0.004 


0.006 


0.007 


0.009 




0.008 


0.007 


0.007 


0.007 


0.006 


0.005 


0.003 


B{D**lu) 


0.022 


0.014 


0.006 


0.000 


0.000 


0.008 


0.006 


B((D^ir) non _ res £v) 


0.024 


0.017 


0.007 


0.004 


0.004 


0.004 


0.004 


D^iv form factors 


0.013 


0.013 


0.012 


0.011 


0.010 


0.008 


0.006 


D**£u form factors 


0.003 


0.002 


0.002 


0.001 


0.001 


0.001 


0.004 


unfolding 


0.015 


0.015 


0.015 


0.015 


0.015 


0.015 


0.015 


binning 


0.001 


0.001 


0.001 


0.001 


0.001 


0.000 


0.001 


efficiency 


0.008 


0.011 


0.012 


0.009 


0.008 


0.005 


0.004 


total 


0.052 


0.041 


0.029 


0.024 


0.022 


0.022 


0.021 



B. Systematic Uncertainties 

The different contributions to the systematic error are shown in Tables IlllffVl The total 
systematic error in Table [III corresponds to the quadratic sum of these components. 

The uncertainties related to the different background components in M\ are estimated 
by varying the respective background normalization factors within ±1 standard deviation. 

We consider both variations of the B — ► D^*'iv branching fractions and form factor 
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TABLE IV: Same as Table [EH for ((M| - (Mj 2 ,)) 2 ). 







A((Ml - 


/ i\/f2 \\2\ 
\ M XI) 1 


[KjG V / 


c ) 




TP* (CnM^ 


U. i 


0.9 


1.1 


1 Q 


1.0 


1 1 
1. 1 


1 o 

i.y 


secondary /fake leptons 


0.167 


0.109 


0.050 


0.023 


0.009 


0.005 


0.002 


combinatorial background 


0.028 


0.018 


0.009 


0.005 


0.003 


0.002 


0.001 


continuum 


0.000 


0.000 


0.000 


0.000 


0.000 


0.001 


0.000 


B — > X u tv background 


U.UU4 


0.004 


0.004 


U.UUo 


n nnQ 
U.UUo 


n nn^ 
U.UUo 


n nn^ 
U.UUO 




0.013 


0.010 


0.007 


0.004 


0.002 


0.002 


0.003 


B(D £u) 


0.216 


0.169 


0.102 


0.049 


0.042 


0.011 


0.009 


vv " ynon — res.^^/ 


0.168 


0.125 


0.058 


0.041 


0.024 


0.004 


0.004 


D^lv form factors 


0.029 


0.028 


0.024 


0.019 


0.017 


0.016 


0.007 


D**tv form factors 


0.013 


0.009 


0.006 


0.003 


0.004 


0.001 


0.004 


unfolding 


0.035 


0.035 


0.035 


0.035 


0.035 


0.035 


0.035 


binning 


0.001 


0.001 


0.000 


0.001 


0.001 


0.001 


0.002 


efficiency 


0.025 


0.032 


0.027 


0.014 


0.013 


0.005 


0.002 


total 


0.327 


0.244 


0.137 


0.080 


0.064 


0.040 


0.038 



TABLE V: Same as Table HI] for (Af£). 



A(i\4) (GeV 4 /c 8 ) 



^min (GeV) 


0.7 


0.9 


1.1 


1.3 


1.5 


1.7 


1.9 


secondary /fake leptons 


0.46 


0.31 


0.16 


0.09 


0.04 


0.02 


0.00 


combinatorial background 


0.08 


0.05 


0.04 


0.03 


0.01 


0.01 


0.00 


continuum 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


B — > X u £u background 


0.04 


0.04 


0.04 


0.04 


0.05 


0.06 


0.07 




0.07 


0.07 


0.06 


0.05 


0.04 


0.03 


0.02 


B(D**£u) 


0.41 


0.30 


0.14 


0.04 


0.04 


0.05 


0.06 


B(( J DW 7 r) non _ res .^) 


0.38 


0.28 


0.12 


0.08 


0.06 


0.04 


0.04 


D^lu form factors 


0.15 


0.14 


0.13 


0.12 


0.10 


0.09 


0.04 


D**lv form factors 


0.04 


0.03 


0.02 


0.01 


0.01 


0.00 


0.03 


unfolding 


0.16 


0.16 


0.16 


0.16 


0.16 


0.16 


0.16 


binning 


0.01 


0.01 


0.01 


0.00 


0.00 


0.01 


0.00 


efficiency 


0.09 


0.13 


0.13 


0.09 


0.08 


0.04 


0.03 


total 


0.77 


0.58 


0.36 


0.26 


0.23 


0.21 


0.20 
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FIG. 3: Graphical representation of the results in Table HH The error bars indicate the statistical 
and total experimental errors. 



shapes. For the former, the ranges of variation are taken from Ref. 3l|. For the latter, the 
curvature p 2 in the form factor parametrization 16j is varied within 1.56±0.14 (1.15±0.16) 
for B — > D*lv (B —>■ Dlv) [32]. For B —>■ D*iv^ we also vary the form factor ratios R\ and 
R 2 M 



The LLSW model 18[ predicts the relative abundance and the form factor shape of the 
different components in B — > D**lv only. To obtain the absolute branching fractions of 
the B -»• D**£u components and of B -> (L'(* ) 7r) non _ rcs /z/, we use -> Z^+z/) = 

(5.6 ± 1.6) x 1(T 3 [3lj, the recent Belle measurement of 8{B -> flWf&j and the total 
semileptonic branching fraction 3l|. The uncertainty assigned to the B — > D**iv branching 
fractions in Tables IHIfiVl reflects the uncertainty in these measurements and the change in 
the -B — > D**£u composition when varying the LLSW parameters within their allowed range. 
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TABLE VI: Correlation coefficients between (Ml) measurements. 



ZT 1 * 
min 

(GeV) 


0.7 


0.9 


1.1 


{Ml) 
1.3 


1.5 


1.7 


1.9 


0.7 


1.000 


0.932 


0.786 


0.615 


0.481 


0.168 


0.071 


0.9 




1.000 


0.888 


0.715 


0.573 


0.241 


0.116 


1.1 






1.000 


0.849 


0.693 


0.363 


0.194 


(M 2 X ) 1.3 








1.000 


0.804 


0.470 


0.254 


1.5 










1.000 


0.591 


0.308 


1.7 












1.000 


0.363 


1.9 














1.000 



TABLE VII: Correlation coefficients between (Af|) and {(Ml — {Ml)) 2 ) measurements. 



7?* 

min 






{{Ml 


- {Ml)) 2 ) 






(GeV) 


0.7 


0.9 


1.1 


1.3 


1.5 


1.7 


1.9 


0.7 


0.897 


0.847 


0.788 


0.713 


0.576 


0.306 


0.102 


0.9 


0.777 


0.843 


0.804 


0.726 


0.608 


0.356 


0.144 


1.1 


0.548 


0.615 


0.757 


0.690 


0.606 


0.426 


0.211 


{Ml) 1.3 


0.328 


0.371 


0.483 


0.718 


0.599 


0.476 


0.260 


1.5 


0.223 


0.263 


0.346 


0.481 


0.702 


0.559 


0.280 


1.7 


-.051 


-.031 


0.035 


0.126 


0.237 


0.846 


0.296 


1.9 


-.060 


-.047 


-.007 


0.040 


0.075 


0.228 


0.865 



The SVD algorithm used to unfold the measured M| distribution requires the detector 
response matrix, i.e., the distribution of measured versus true values of M\. We determine 
this matrix from the MC simulation. To study the systematics related to unfolding and a 
possible mismodeling of the detector response, we change the amount of bin-to-bin migration 
by varying the effective rank of the detector response matrix, the main tunable parameter 
of the SVD algorithm. We have further studied a change of the binning of the unfolded 
distribution and the effect of disabling the bin-to-bin efficiency correction. 

C. Correlations 

Due to overlapping events, the moment measurements corresponding to different lepton 
energy thresholds are highly correlated. Systematic uncertainties are another source of 
correlation. We have estimated the correlations due to both sources using a toy MC approach 
based on 50,000 simulated measurements. The results for the self- and cross-correlation 
coefficients are given in Tables IVIflXl 
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TABLE VIII: Correlation coefficients between (Ml) and (Ml) measurements. 



E* ■ 

mm 








(Ml) 








(GeV) 


0.7 


0.9 


1.1 


1.3 


1.5 


1.7 


1.9 


0.7 


0.983 


0.933 


0.830 


0.683 


0.523 


0.194 


0.073 


0.9 


0.890 


0.974 


0.910 


0.765 


0.606 


0.264 


0.117 


1.1 


0.704 


0.810 


0.976 


0.857 


0.707 


0.380 


0.196 


{Ml) 1.3 


0.508 


0.601 


0.774 


0.980 


0.800 


0.479 


0.258 


1.5 


0.383 


0.469 


0.614 


0.764 


0.985 


0.597 


0.305 


1.7 


0.079 


0.137 


0.273 


0.403 


0.539 


0.994 


0.357 


1.9 


0.017 


0.052 


0.136 


0.208 


0.271 


0.348 


0.995 



TABLE IX: Correlation coefficients between ((M| - (Ml)) 2 ) measurements. 



E* ■ 

mm 








mi 


-(Ml 


» 2 > 






(GeV) 




0.7 


0.9 


i.i 


1.3 


1.5 


1.7 


1.9 




0.7 


1.000 


0.939 


0.838 


0.698 


0.534 


0.167 


-.024 




0.9 




1.000 


0.901 


0.732 


0.586 


0.195 


-.011 


mi- 


1.1 






1.000 


0.793 


0.638 


0.262 


0.034 


-(Ml)) 2 ) 


1.3 








1.000 


0.731 


0.340 


0.102 




1.5 










1.000 


0.484 


0.146 




1.7 












1.000 


0.296 




1.9 














1.000 



TABLE X: Correlation coefficients between (Ml) measurements. 



E* ■ 

mm 

(GeV) 


0.7 


0.9 


1.1 


(Ml) 
1.3 


1.5 


1.7 


1.9 


0.7 


1.000 


0.932 


0.784 


0.601 


0.442 


0.111 


0.017 


0.9 




1.000 


0.877 


0.684 


0.524 


0.168 


0.051 


1.1 






1.000 


0.817 


0.651 


0.297 


0.137 


(Ml) 1.3 








1.000 


0.780 


0.421 


0.212 


1.5 










1.000 


0.557 


0.270 


1.7 












1.000 


0.346 


1.9 














1.000 
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IV. SUMMARY 



We have measured the first, (M|), and the second central and non-central moments, 
((M x — (M x )) 2 ) and (M x ), of the hadronic mass squared spectrum in B — > X c lv decays for 
lepton energy thresholds ranging from 0.7 to 1.9 GeV. Using a toy MC approach, we have 
also evaluated the full covariance matrix for this set of measurements. 

It is expected that this measurement, combined with measurements of the semileptonic 
branching fraction, moments of the lepton energy spectrum in B — > X c iv decays and possibly 
other moments, will lead to an improved determination of 6-quark mass m;, and the CKM 
matrix element \V c b\ @, 0, 0, [5| . 
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